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Abstract 

We examine applications of polynomial Lie algebras slpd{2) to solve physical tasks 
in Gj„^-invariant models of coupled subsystems in quantum physics. A general operator 
formalism is given to solve spectral problems using expansions of generalized coherent 
states, eigenfunctions and other physically important quantities by power series in 
the slpd{2) coset generators V±. We also discuss some mappings and approximations 
related to the familiar sl{2) algebra formalism. On this way a new closed analytical 
expression is found for energy spectra which coincides with exact solutions in certain 
cases and, in general, manifests an availability of incommensurable eigenfrequencies 
related to a nearly chaotic dynamics of systems under study. 

1 Introduction. General remarks 

Recently, in a series of our papers /1-4/ a new efficient Lie- algebraic approach has been 
suggested to solve both spectral and evolution problems for some nonlinear Gmti-invariant 
models of coupled subsystems in quantum physics. It was based on exploiting a formalism 
of polynomial Lie algebras Qpd as dynamic symmetry algebras g^^ of models under study, 
and besides generators of these algebras gpd can be interpreted as Gm^j-invariant "essen- 
tial collective dynamic variables" in whose terms model dynamics are described completely. 
Specifically, this approach enabled to develop some efficient techniques for solving physical 
tasks in the case g^^ = slpd{2) when model Hamiltonians H are expressed as follows 

H = aVo + gV+ + g*V. + C = AY + C, [Va,C]=0, V. = (LI) 

where C is a function of Gj„^-dependent model integrals of motion Ri and Vq, V± are the 
slpd{2) generators satisfying the commutation relations 

[Vo, v±] = ±r±, [v., V+] = MVo + 1) - MVo), (1.2) 

with the structure function i'niyo) being a polynomial ipniVo) = ^{{Rj}) IYi=i(yo~ ^i{{Rj})) 
of the degree n in Vq. 

For example, the three-boson model Hamiltonian 

II2 = Uiafai + u;2«2''^2 + ^3(^3^3 + fl'('^^'^2^)'^3 + 9*{cii(^2)ci3 (1-3) 
can be expressed in the form (1.1) if using the substitutions 

Vq = {Ni + N2- A^3)/3, V+ = (a]^a^)a3, a = ui + uj2 - UJ3, Ni = afai, 
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2C = i?i(a;i - 002) + R2{uJi +UJ2 + 2uj^), = Ni - N2, 3R2 = Ni + N2 + 2N3 (1.4a) 
In this case the structure function '0„(Vo) = ip3{Vo) is given as follows 

MVo) = \i2Vo + R2- i?i)(2Vo + Ri + R2)i-Vo + R2 + I) (1.46) 

All techniques developed were based on using expansions of physically important quan- 
tities (evolution operators, generalized coherent states (GCS), cigcnfunctions etc.) by power 
series in the slpd{2) coset generators V±. Besides, in the Schroedinger picture one has ex- 
ploited decompositions 

L{H) = J2mi]): {V+V^-MVo))\Lim)^0 (1.5) 

e 

of Hilbert spaces L{H) of quantum states of model in direct sums of the subspaces ([/•«]) 
which are irreducible with respect to joint actions of algebras slpd{2) and groups G'lnv and 
describe specific "s/p(i(2)-domains" evolving independently in time under action of the Hamil- 
tonians (1.1). The subspaces L{\1^) are spanned by basis vectors 

\[%V) = [{Mh-rv))^^^\-"'Vl\[k\), (V'n(x))(^) = nVn(x-r), 

Vo\[%v) = {h + v)\[%v), Ri\U,'^)-k\U,v), i^n{Ro)=Myo)-V+V^ (1.6) 

where is the lowest vector (\^_|[/j]) = 0, ■0„(Zo) = 0) of L{\li\). 

Then, using Eqs. (1.1)-(1.2) one may get Jacobi-type three-term recurrence relations for 
amplitudes Qy{Ef) = {[lj\\v\Ef) of expansions of energy cigcnstates \Ef) in bases 
Besides, energy spectra {-E/} of bound states arc given by roots of certain spectral functions 
(polynomials for the compact version slp^i^) = sUpd{2)) which are determined for given 
structure functions ipni^) with the help of similar recurrence relations /1-3/. Another way, 
exploiting the Bargmann-type representation of the slpd{2) generators, 

V-^ = z, Vq = zd/dz + lo, V_ = z~^il!n{zd/dz + /q), (1.7a) 

reduces these tasks to solving some singular differential equations /1-3/. When using a 
conjugate to (1.7a) representation of the slpd{2) generators, 

V_ = d/dz, Vq = zd/dz + Iq, V+ = i/j^izd/dz + lo){d/dz)~'^ (1.76) 

this way leads to solving the Riccati-type equations for structure functions ip3{x) of the 
degree n = 3 (that is the case for the Hamiltonian (1.3)). In the paper /4/ some integral 
expressions were found for amplitudes Qy{E), eigenenergies {Ea} and evolution operators 
Unit) with the help of a specific "dressing" (mapping) of solutions of some auxihary exactly 
solvable tasks with the dynamic algebra s/(2). 

However, all these and other results do not yield simple working formulas for analysis 
of models (1.1) and revealing different physical effects (e.g., collapses and revivals of the 
Rabi oscillations /5/) at arbitrary initial quantum states of models. Besides, solutions /4/ 
of spectral tasks manifest so-called "quantum discontinuities" /6/: a disappearence of wave 
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functions when attaining the hmit of auxihary s/(2) Hamiltonians that makes difficult to 
compare completely quantum models with their semi-classical analogs. Therefore, it is nec- 
essary to develop some simple techniques enabling to display important physical peculiriaties 
of models (1.1)-(1.2). In the case g^^ = sl{2), when the structure functions il^n{x) = 'ip2{x) 
are quadratic functions ip2{x) = {j =t x){^^j + 1 — x), Iq = =Fj)) the GCS formalism of the 
group orbit type /7/ is known to be an efficient tool for analyzing both linear /7/ and non- 
linear /7-9/ models. This formalism based on properties of the SL{2) group displacement 
operators Sv{C) — 6xp(^y_)_ — i*VJ), ^ — rexp{—i9) yields exact solutions /7/ for linear 
models and variational schemes (corresponding to the Ehrenfest theorem) to obtain effective 
mean-field approximate solutions for non-linear models /8-9/. 

Below we examine some possibilities of generalizations of this formalism for solving spec- 
tral problems of models (1.1)-(1.2) (Section 2) and give a variational scheme to find "smooth" 
sZ(2)-approximations of these solutions (Section 3) using an isomorphism of the slpd{2) alge- 
bras to special subalgebras of the extended enveloping algebra U (s/(2)) of the familiar algebra 
s/(2). This isomorphism is established via a generalized Holstein-Primakoff mapping given 
as follows /1-3/ 



Yo = Vo-loTJ, J = L Y+ = V+ 



{jTYo){±j + l + Yo) _ 
2' - MVo + 1) ' - " ^ ' ^ 

where Ya are the sl{2) generators, upper and lower signs correspond to the su{2) and su{l, 1) 
algebras respectively. In Section 4 some prospects of further studies along these lines are 
briefly outlined. 



2 A general operator formalism to solve spectral prob- 
lems 

As is known /7/, the Hamiltonians (1.1) are simply diagonalized with the help of operators 

SviO = eMCV+ - CV-) = exp[i(r)e^V+] exp[-2 In c{r)Vo\ exp[-i(r)e-^V_], ^ = re'' 

(2.1) 

when Va arc generators of the familiar s/(2) algebra (t(r) = tanr, c(r) = cosr for su{2) and 
t{r) = tanh r,c(r) = coshr for su{l, 1)). Indeed, using the well-known sl{2) transformation 
properties of operators Va one finds the transformation 

H H{0 = SvimSviO^ = C + VoAoia, g; + V+A+{a, g; + V.Al{a, g; (2.2a) 

of the Hamiltonians (1.1) under the action of operators iS'v'(.^). Then, supposing A^{a, g; ^) = 
we find a value of the parameter ^ diagonalizing the Hamiltonian H[^). For example, 
in the case of the su{2) algebra we have /7/ 

H{^o) = S{^o)HS{^oy + Vo^o? + A\g\\ = ^arci^^, (2.2) 

z\g\ a 

and the corresponding eigenenergies E{\1^]V]^q) and eigenfunctions |[/i];'y;Co) are expressed 
as follows 

a) E{[k]-v-^o) ^C+{-3+v)^a^ + A\g\\ (2.3a) 
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b) m;v;^o) = Sv{^oy\[h];v;)=exp{-^oV+ + CoV-)m;v-) = (cosV^^x 
E -v + 2j + l-J-v + l; sin^ /; ) (2-36) 

where F{...) is the Gauss hypergeometric function. An equivalent way /8/ to obtain the 
results (2.3) is based on using the stationarity conditions 

dE([k];v;0 dE([k];v;0 _^ 

89 ' dr 

for the energy functional E{[lj\]v; ^) — {\li\;v] ^\H\[l^]v; ^) defined with the help of the SL{2) 
GCS = v; ) as trial functions. 

Both ways above essentially exploit the finite-dimensionality of the s/(2) adjoint (vec- 
tor) representation (cf. Eq. (2.1a)) and well-known (due to Eq. (2.1)) explicit expansions 
of the SL{2) GCS in orthonormahzed basis states. However, for polynomial Lie algebras 
slpd{2) the situation is more complicated since their adjoint representations defined by re- 
peated commutations of arbitrary elements arc infinite-dimensional as it follows from Eq. 
(1.2). Furhtermore, GCS exponential operators SyiO — cxp(.^K|_ — ^*V"_) have not explicit 
expressions for matrix elements in orthonormahzed bases (1.6) as these exponentials are not 
elements of Lie groups but only correspond to quasigroups (pseudogroups) /lO/ which have 
no simple analogs of the "disentangling theorem" (2.1) providing expansions of operators 
SviO finite products of one-parameter subgroups /lO, 11/. Therefore, in this case a 
direct generalization of results (2.3) is impossible. 

Nevertheless, taking into account Eqs. (1-2), (1.8) one may apply the diagonalizing 
scheme (2.2) using represntations of diagonalizing operators S{^) by power series 

oo 

S{0- yISf{Vo;0, V-' = V^{[MVo)f^)-\k>0 (2.5) 

/=-oo 

with undetermined (unhke those for the sl{2) algebra - cf. (2.1) and (2.3)) coefficients 
5'/(Vo; ^) (which, when being known, provide possibilities of explicit calculations of any physi- 
cal quantities with the help of Eqs. (1.2), (1.6)). For diagonalizing operators S{^) = Sv{C) — 
exp{^V-^ — C,*V-) (if they exist) these coefficients may be taken in the form Sf{Vo;$,) — 
exp(i/6')(T/(Vo; r), ^ = rexp(i6'), and satisfy the equations 

^^^^^ - aj.^iVo; r) + ^„(yo + /)a/+i(l^o; r) = 5(r)5;,o, / = 0, 1, ... (2.6) 

whose solutions may be represented by power series in r (via direct expansions of exponents 
Sv{C)) or obtained in an integral form with the help of the "sZ(2) dressing" procedure /4/. 
In general cases these coefficients satisfy the equations 

oo 

E [^(Vo)]^^^WK)-/;e)^K^o-/; = 4,0 (2.7) 

/=-oo 

following from the unitarity conditions SS"^ = S'^S = I. 
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Then, substituting Eq. (2.5) in the scheme (2.2) one gets after some algebra nonhnear 
analogs of Eqs. (2.2) 

oo 
f=-oo 

oo 

hfiVo;0= E mVo)f-^^S,{Vo + f-k;OHVo + f-k)S:_f{Vo + f-k;0+ 

fc=— oo 

gi;{Vo + / - k)Sl^,^fiVo + / - - 1; + /5fc-i-/(K, + / - A: + 1; 0], 

h.f{Vo;0 = mVo-f;OmVo)]^^\ mVo)t^^^mVo + f)]^^Y', />0, (2.8a) 

b) H{^o) = S{^o)HS{^oV ^C + ho{Vo;^o): E{[k];v;0 ^ C + {[k]; v\ho{Vo; OM; v) (2.86) 

expressed in terms of the coefficients Sf{Vo;^) (hereafter the subscript n in '0n(H) will be 
omitted for the sake of the notation simplicity). As is seen from Eq. (2.8b) the diagonalized 
Hamiltonian i?(Co) has (unlike (2.2b)) an essentially non-linear dependence in Vq determined 
by coefficients Sf{Vo;^o) which satisfy (additionally to Eqs. (2.7)) the operator recurrence 
relations following from the condition S{^o)H = H{^o)S{^o)), 

SfiVo; Co)[aVo - ho{Vo + /; 6)] + 9Sf-i{Vo + 1; ^o) + 9*Sf+i{Vo - 1; Co) = 0, / = 0, ±1, ±2, .., 

oo 
n=—oo 

gi;{Vo - n)S:^,iVo - n - 1; Co) + 9*S:_,{Vo -n + 1; Co)] (2.9a) 
or the operator equations 

oo 

0= E mVo)]^''-f^Sn{Vo + f-n;^oMVo + f-n)S:_f{Vo + f-n;^o) + 

n=— oo 

g^{Vo + f-n)S:^,_f{Vo + f-n-l;Co)+9*S:_,_fiVo + f-n + l;Co)], / = ±1,±2,.. (2.96) 

resulting from the condition /i/(Vo;Co) = 0, / = ±l,ih2,.. (a direct generalization of the 
condition A_^{a,g;^) = in (2.2)). Note that Eqs. (2.9), in general, determine both a 
suitable functional form of Sf{Vo;^) and a value of the parameter ^ diagonalizing the 
Hamiltonian H(^)- 

So, the formalism of the slpd{2) algebras enabled to to get a general operator scheme of 
diagonalizing the Hamiltonians (1.1) with the help of solving the (infinite) set of algebraic 
operator equations (2.7)-(2.9). Evidently, without using some specifications of diagonalizing 
operators S{^) the task of solving these equations is equivalent to that for finding amplitudes 
Q,{Ef) = {[h];v\Smh];f) = [(^„(/o + /))(^V(^n(/o + t^))(^)]'/'<5)_,(/o + t'; Co) as Eqs. (2.9b) 
resemble those for Q^{Ef). Note that in the case of the compact sUpd{2) algebra, having 
finite-dimensional (with dimensions equal to (i([/i])) irreducible subspaces L{[li]), it is possible 
to simplify the task restricting oneself by the consideration of Eqs. (2.7)-(2.9) on each 
L([Zj]) independently. Then, due to the relation (V±)''(['*l^+^ !£,([;.]) = all series in Eqs. 
(2.5), (2.7)-(2.9) are terminating. Specifically, wave eigenfunctions \Ef) — S^lk]; f) may 
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be represented by polynomials \Ef) = A{Vo)Y[j{V+ — ■^i(^)) and the energy functional 
E{[h]'-,f)s = {[h]', f\SHS^\[li]; f) are written down in the form of a sum of d{[li]) spectral 
functions as it is prescribed for such classes of models by the algebraic Bethe ansatz /1 2/. 
(In essence, we obtain in such a manner a new formulation of this ansatz for a wide class 
of models in terms of the sUpd{2) algebras which is simpler and more efficient (cf. /3/) in 
comparison with its initial version /12/.) 

However, even such simphfications enable us to get simple closed expressions only for 
little dimensions d{[li\) of the sUpii{2) irreducible subspaces -^([^j])- At the same time many 
physical quantum states of models (1.1), e.g., such as coherent and squeezed states in models 
(1.3), have non-zero projections on all subspaces L([/j]). Therefore, for physical applications 
it is necessary to get some closed expressions like Eqs. (2.3) for energy eigenvalues and 
wave eigenfunctions which would describe main features of model dynamics with a good 
accuracy. One exampl of such analytical approximations was obtained in /1-3/ by mapping 
(with the help of the change ^ Hamiltonians (1.1) in Hamiltonians 5.1(2) which 
are linear in s/(2) generators and have on each fixed subspacc i^([/i]) equidistant energy 
spectra given by formulas like Eq. (2.3a) (but with modified constants a,g). However, this 
(quasi)equidistant approximation is suitable for little or very big dimensions d{[li]) and does 
not enable to display many peculiarities (e.g., availability and a fine structure of collapses and 
revivals of the Rabi oscillations) of models (1.1). Therefore, below we describe an alternative 
approximation applying the variational scheme (2.4) with SL{2) GCS as trial functions to 
Hamiltonians (1.1) expressed with the help of Eqs. (1.8) as functions of sl{2) generators Fq. 

3 A variational scheme of determining energy spectra 
with the help of 6'L(2)-coherent states 

Hamiltonians (1.1) re- written in terms of Ya have the form 



H^aYo + gY+ 



MVo + 1) 

+ 9 



+ -Y_ + C + a{±j + k) 



'\ ij T Yo){±j + l + Yo)^' \{JT Yo){±j + I + Y0) 

. . ^^-^'^ 

which is essentially non-linear in sl{2) generators Y^. Therefore, in general, it is unlikely to 
diagonalize them with the help of operators Sy{^) — exp(^l+ — ^*y_). However, it is natural 
to apply associated with these operators SL{2) GCS 

\Uv;O-SY0m;v;)^eM-^y+ + CY-)m;v;), (3.1) 

as trial functions in the variational scheme (2.4) that results in non-linear analogs of Eq. 
(2.3a) for approximate energy eigenvalues. Such an approximation may be called as a 
"smooth" sl{2) approximation since it, in fact, corresponds to picking out a "smooth" (due 
to analytical nature of SL{2) group elements) s/(2) factor exp(^y_|_ — $,*Y-) in the exact 
diagonalizing operators S{^). 

Specifically, application of this procedure to Hamiltonians with the sUpd{2) dynamic 
symmetry yields after some algebra the following expressions 

E{[%v, ^o) = C + a{k + j) + a{-j + v) cos2r - 2\g\{cos^ rfU-v){^^l^_^ ^ 



7)1 

f>0 
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Et(k,rJ) 



(t^r)2(/--)+i(/ + l)! 



+ 1 + /) 



X 



/-i)!^(2j-/)(/ + i; 

f 1; sin^ r)F{-v, -v + 2j + 1; f - v + 2; sin^ r) (3.2) 
re*^) where e*^ = due to the second condition (2.4) 



(f - v)\if + 1 - v)\{2j 

F{-v,-v + 2j + l-J-v 

for energy eigenvalues E{[li]; v; = 
and diagonahzing values of the parameter r are determined from solving the equations 







.2/ 



- |4« J - (l + a )(2/+ 1)] J p-^^^ 



}, a = -tgr (3.3) 



resulting from the first condition (2.4). 

As is seen from Eq. (3.2), spectral functions E^Qq, j; f) are non-linear in the discrete 
variable v labeling energy levels that provides a non-equdistant character of energy spectra 
within fixed subspaces -^^([/i]) at d{[li]) > 3; besides, due to the availability of square roots in 
expressions for these functions different eigenfrequencies ujy = E^/h are incommensurable: 
mujy^ 7^ na;„2 that is an indicator of an origin of collapses and revivals of the Rabi oscillations 
/5/ as well as of prechaotic dynamics. We also note that Eqs. (3.2)- (3.3) give exact results 
at little dimensions d{[li]). 



4 Conclusion 

So, we have obtained a general operator scheme for diagonahzing Hamiltonians (1.1) and a 
smooth approximation for solutions of its defining equations with the help of the mapping 
(1.8) and the variational scheme (2.4) using the SL(2) GCS as trial functions. This approx- 
imation may be used as an initial one in iterative schemes of solving Eqs. (2.9) (re-written 
in "the si {2) languauge") which are similiar to those developed to examine non-linear prob- 
lems of classical mechanics and optics /II/. Further investigations along this hne may be 
also related to a search of suitable specifications of the operators S{^) (besides the form 
S{$,) — Sv{C,)) reducing solutions of Eqs. (2.9) to determining a value providing an ex- 
act or a sufficiently accurate approximation for diagonalization of the Hamiltonian (1.1) in 
scheme (2.3a). Another way to develop the results above concerns some simplifications of 
the formulas (3.2) via using different properties, including integral representations, of the 
hypergeometric functions F{a,b;c;x). For the case of ■03(^) it is also of interest to com- 
pare results of such approximations with exact calculations obtained by considering exactly 
solvable cases of the Riccati equations yielded by the slpd{2) representation (1.7b). 

Finally, general ideas of the analysis above may be extended to solve evolution problems. 
Specifically, a version of general operator formalism for these tasks was formulated in /4/, 
and a version of obtaining a variational dynamics in the mean-field approximation can be 
found following the approach of the paper /8/. Namely, looking for the evolution operator 
Unit) (with H given by Eq. (1.1')) in the form Unit) = exp{—z{t)Y+ + z{t)*Y^) and using 
the sl{2) GCS \z{t)) = cxp(— 2;(t)y+ + z(t)*Y^)\iJjo) as trial functions in the time-dependent 
Hartree-Fock varitional scheme with the Lagrangian £ = {z(t)\{id / dt — H)\z{t)) one gets 
the sl{2) Euler-Lagrange equations 

. m . m 
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for "motion" of the s/(2) GCS parameters; here Ti = {z{t)\H\z{t)) and p = j cos6, q = (f),z = 
tan(6'/2) exp(— i0) for sm(2) and p = j cosh. 9, q = (f),z = ta.nh{9/2)exp{—i(f)) for sm(2). 
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